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ABSTRACT 


In this thesis, we emulate a Wind Turbine Generator by driving a Doubly Fed Induction 
Generator (DFIG) via a DC motor with variable input torque capability. The two circuits 
of concern are the DFIG and Supply-side circuits. They are electrically coupled with 
back-to-back Space Vector Modulators (SVMs) that are coupled through a DC Link 
consisting of a capacitor bank. The SVMs, with DC Link, provide bi-directional power 
flow between the DFIG rotor and Supply-side power supply. 

The Supply-side circuit senses voltage, current and DC Link voltage (VDC), and 
sends the sensed inputs to a Field Programmable Gate Array (FPGA), which is used to 
derive control of VDC so that it is maintained at 200V via SVM. 

The DFIG circuit senses rotor current, stator voltage, rotor speed, and rotor 
electrical position, which are used to derive a current reference that is used to control the 
rotor speed using a Proportional Integral Controller (PI). For our application we show 
that the DFIG rotor speed can be controlled (subsynchronous, synchronous, or 
supersynchronous) for a given input torque and thus distributes the power generated by 


the DFIG as predicted between the rotor and stator. 
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EXECUTIVE SUMMARY 


In this thesis, we emulate a Wind Turbine Generator by driving a Doubly Fed 
Induction Generator (DFIG) via a DC motor with variable input torque capability. 
Previous research has shown we can use a DFIG with back-to-back SVMs to accomplish 
a decoupled Supply-side, and DFIG-side, control scheme while allowing power flow to 


occur between the two systems. 


The two circuits of concern are the DFIG and Supply-side circuits. The DFIG and 
Supply-side circuits are electrically coupled with back-to-back (DFIG-side and Supply- 
side) Space Vector Modulators (SVMs) that are coupled through a DC Link consisting of 
a capacitor bank. The back-to-back SVMs, with DC Link, provide bi-directional power 
flow between the DFIG rotor and Supply-side power supply. Bi-directional power flow is 
achieved when the DFIG is controlled at supersynchronous speeds, and input torque is 
enough to overcome losses, and subsynchronous speeds (Supply-side provides power to 


the rotor). 


The Supply-side power supply (3-phase AC) provides power to the Supply-side 
circuit (stepped down by a VARIAC) and the Stator of the DFIG. The Supply-side circuit 
senses the supply side voltage, current and DC Link voltage (VDC) and sends the sensed 
inputs to a Field Programmable Gate Array (FPGA). The FPGA is programmed via a 
Simulation software package called Simulink to achieve the desired control of the 
specified parameters. For the Supply-side circuit we use the mentioned sensed inputs to 
derive control of VDC so that it is maintained at 200V via SVM. There is also a built-in 
protection feature in the Supply-side circuit that will turn off the SVM if a VDC of 240V 


is reached. 


The DFIG circuit power supply is provided by the DC Link and DFIG-side SVM. 
The DFIG circuit senses rotor current, stator voltage, rotor speed, and rotor electrical 
position (as provided by an encoder). These inputs, like the Supply-side circuit, are then 
sent to an FPGA which is programmed to control the DFIG parameters. The DFIG circuit 


will take the rotor speed input and derive a current reference that is used to control the 


XV 


rotor speed using a Proportional Integral Controller (PI). It is important to show that we 
can control the DFIG rotor speed since traditional Wind Turbine Generators would use 
turbine blade pitch control in tandem with a synchronous generator to provide a constant 
speed, and hence a constant output frequency. Further, it is noteworthy that a wind-speed 
profile can be adapted to a specific Wind Turbine (outfitted with a DFIG) in order to 
program the DFIG to run at specified speeds for a given wind speed. This tool can be 
used to provide an optimized rotor speed for a given wind speed. For our application we 
will show that the DFIG rotor speed can be controlled (subsynchronous, synchronous or 
supersynchronous) for any given input torque and thus distributes the power generated by 
the DFIG as predicted between the rotor and stator. Further, this experiment will 
demonstrate the bi-directional power flow of the back-to-back SVMs showing power 


recovery via the DFIG rotor. 


As illustrated in the below figure, the Supply-side SVM converts a 3-phase AC 
source to a constant DC Link voltage of 200V. The DC Link voltage is then inverted into 
a 3-phase waveform via the DFIG-Side SVM (during modes for which power is being 
delivered to the DFIG). The digital control for the Insulated Gate Bipolar Transistor 
(IGBT) network is implemented via the FPGA. The FPGA interface samples the rotor 
current and the stator voltage of the DFIG through an analog to digital (A/D) converter. 
There is an input from the encoder, 9,, to provide the position and speed of the rotor to 
the FPGA through the interface circuit board. An algorithm then calculates the slip 
frequency of the DFIG which 1s used to derive the needed rotor current frequency and 
amplitude to maintain a specified rotor speed, and therefore a controllable output 


frequency. 
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Simplified Circuit of Wind Turbine DFIG System. 


The below figure demonstrates an input torque transient representing a large wind 
gust. By commanding a specified supersynchronous speed, while taking into account 
windage and friction losses of the system, it is clear that the input torque transient 
produces a power flow reversal in the rotor. In other words, the rotor mode transitions 
from power flow from the rotor (negative rotor current) to power flow to the rotor 


(positive rotor current). 


In addition, this transient plot (simulation model) shows a smooth transition as the 
power flow reverses between power flow from the rotor to power flow to the rotor. This 
transition 1s modeled using a constant commanded rotor speed along with a 10% decrease 
in input torque. The next two figures (experimental results) parallels this transition with 
the steady state captures corresponding to the simulation transient. The experimental 
results were run through a low pass filter with a corner frequency of 80 Hz in order to 


create a useful plot of the event. Due to windage and friction losses, a commanded 


XVII 


supersynchronous speed was ordered to achieve the power flow transition results during 
the constant rotor speed experiment. If the experimental results are closely examined, it 1s 


interesting to note the small, low frequency ripple due to speed control beating. 


Translent Gllh constant rotor qpeedanda 10% decrease in Lagut tongue 
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Simulation: Input Torque Transient Showing Rotor Power Reversal. 
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Experimental: Rotor Current and Voltage with Power Flow from the Rotor. 
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Experimental: Rotor Current and Voltage with Power Flow to the Rotor. 


In conclusion, we showed that we can control the DFIG rotor excitation to 
provide a constant output frequency. Our application demonstrated that the DFIG rotor 
speed can be controlled (subsynchronous, synchronous or supersynchronous) for a given 
input torque and thus distributes the power generated by the DFIG as predicted between 
the rotor and stator. Finally, this experiment also demonstrates the bi-directional power 


flow of the back-to-back SVMs showing power recovery via the DFIG rotor. 
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I. INTRODUCTION 


A. MISSION 


In this thesis, we emulate a Wind Turbine Generator by driving a Doubly Fed 
Induction Generator (DFIG) via a DC motor with variable input torque capability. As 
discussed in [1-4] we can use a DFIG with back-to-back SVMs to accomplish a 
decoupled Supply-side, and DFIG-side, control scheme while allowing power flow to 
occur between the two systems. The control schemes developed with this research can 


branch into shipboard application as discussed below. 


The benefits of researching and developing a Wind Turbine Generation System, 
refitted with a DFIG, span a wide range of purposes for the United States as a whole, as 
well as the United States military. In particular, there are several naval warships still 
running large, separately supplied, gas turbine generators to provide the ship with 
electrical power to vital loads. Considering this one example, a properly researched and 
developed DFIG with precision characteristic control can not only save space (in a space- 
scarce environment), but also provide a more efficient means of power distribution when 
mechanically coupled to the main propulsion gas turbine shaft line. This, of course, 
requires a variable speed constant frequency control scheme since the propulsion turbine 
shaft speed is variable according to the desired ship speed [5]. The applications of the 
DFIG are many and complex, but the simple example above shows a common sense 


approach to save space and fuel consumption for our United States military. 


B. OBJECTIVE 


The goal of this thesis was to build a computer simulation of a Wind Turbine 
Generation System, based on mathematical equations, through the use of Mathworks’ 
Simulink software. Once the simulation was complete, the goal shifted to build a fully 
functioning model of a Wind Turbine Generation System. A series of experiments were 
performed to compare actual results with expected/simulation results. Once the series of 


experiments proved the simulation was an accurate depiction of the model, the machinery 


| 


was put to full test as a Wind Turbine Generation System to be used in future lab work 


and research to hone in features that could be discovered. 


C. APPROACH 


The first step is to generate a computer simulation of the system, shown in Figure 
1. There are a number of ways to simulate the DFIG coupled with a Supply-side, via 
back-to-back SVMs, but the one used for this thesis 1s covered in [1]. Specifically, a 
computer model of the DFIG-side circuit based on mathematical equations 1s simulated in 
the rotor reference frame using dqQ transformation. The advantage of using dq0 
transformation 1s most apparent when attempting to decouple DFIGs real and reactive 
power. By doing so, it 1s possible to separately control generator torque (hence speed) 
and reactive power (hence rotor and stator power factor) when the currents controlled are 
in the slip frequency reference frame. The DFIG-side circuit simulation inputs are input 


torque, stator voltage, rotor current, and rotor speed and position. 
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Figure 1. Simplified Circuit of Wind Turbine DFIG System. 


The Supply-side circuit simulation is based in the stator reference frame using dq0O 
transformation. Not unlike the DFIG-side circuit, there is a distinct advantage in using the 
dqO transformation so as to decouple real and reactive power on the Supply-side circuit. 
This, in turn, allows a separate control of DC Link voltage (via a current reference) and 
reactive power (hence Supply-side power factor). The Supply-side circuit simulation 


inputs are Supply-side voltage, DC Link voltage and Supply-side current. 


The physical system consists of a DFIG, back-to-back SVMs, Student Design 
Center (SDC) [6], an encoder, and a DC motor used to simulate changes in wind speed. 
The DFIG is connected to the grid via the stator windings and to the DFIG-side SVM via 
the rotor windings. The SDC is the mechanism by which the control algorithm 1s 


implemented. The SDC measures the applicable system parameters and outputs the 


appropriate control signals to the SVM to control the rotor currents in a fashion that 
allows for a specific mode of operation. For our application we will show that the DFIG 
rotor speed can be controlled (subsynchronous, synchronous or supersynchronous) for a 
given input torque and thus distributes the power generated by the DFIG as predicted 


between the rotor and stator. 


D. THESIS ORGANIZATION 


This thesis is organized with Chapter II covering the theory of induction 
machines; to include the equations that represent voltage, flux linkage, and reference 
frame transformations and how the voltage and flux equations and machine parameters 
are used to model the system. Chapter III is dedicated to the control algorithm. It explains 
how an algorithm is formulated, based on the equations derived in Chapter II, and 1s 
manipulated to control the simulation and the physical system. Chapter IV displays the 
results of the working system, and compares the actual results to the simulation results. 
Chapter IV also includes the testing procedure used to calibrate and implement the 
encoder and its offset, the method used to transform the DFIG-side rotor currents to a slip 
frequency reference frame, and Xilinx implementation methods. Chapter V discusses the 
conclusions and future research objectives for this thesis. The Appendix contains Matlab 
files, data sheets for equipment used, specific schematic layout of the Simulink systems, 


and transformation derivations used in this thesis. 


Hi. DFIG THEORY 


A. INDUCTION MACHINE EQUATIONS 


Simulink is a powerful tool that can be implemented in the simulation of any 
given system that can be mathematically represented. To understand the simulation of the 
Wind Turbine Generation System, refitted with a DFIG, there must be a thorough 
understanding of the equations that represent the DFIG. The equations used to 
characterize the DFIG are obtained from [7]. The induction machine voltage equations 


are expressed as 


te = ly, ee + Pp | Aaah 
Maher = lv, Le + Pp xe 


where the subscript s and r are denoting the stator and rotor associated variables and A 


(1) 


represents the flux linkage. Finally, p is used as the derivative operator. For a 


magnetically linear system, flux linkage can be represented as equation (2). 


Aabes L, L, Labes 
= ih ; (2) 
A aber (L,, ) L. lober 
where L,, represents the mutual inductance between the stator and rotor. The stator, rotor 


and mutual winding inductances are 


| ey ee ey 0 


2 “~ms 2 “~ms 
L, = ace oe L, a oem a oe 
pile 2 Lins L, a Dy 
(3) 
L, + Lis Sy ae ae en 
L, _ als L, ge Ooms oe 
oe ee L, me Ls 
(4) 
cos 0. cos(6.+2£) cos(@, —22) 
is oe cos(0 — 34) cos 0. cos (6, + 22) 
cos(6,+2£) cos(6. —24) cos 8. 
(5) 


where the variables L}; and Lms are the leakage and magnetizing inductances of the stator 
windings; conversely, Li, and Lm, are the leakage and magnetizing inductances for the 


rotor windings. Lastly, 0,1s the angular position of the rotor. 


When creating an equivalent circuit, a standard convention is to refer the rotor 


variables to the stator windings using the appropriate turn ratio N. 





(6) 
i ae 
V chen = N. Viney 
(7) 
NG 
A aber i N. A aber 
(8) 
2 
posh es 
r= r ; 
7 
(9) 


This turns ratio referral also applies to the magnetizing and mutual inductances and are 











related by 
Lins = a L,, 
N, 
(10) 
N 2 
EAP 
Ne 
(11) 
N 2 
L, -[ : L, 
N, 
(12) 


N N 
where by defining L’ =—L, and L’,. =| — 
N N 





r r 


2 
L,; the rotor inductance variables can be 


referred to the stator windings as shown 1n equations (13) and (14). 
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cos 8. cos(6,+2£) cos(@, —2£) 


L'. =L,, | cos(0, —2£) cos @. cos(6, +24) 


3. 


cos (0, + 22.) cos (0, — 32) cos 0. 


3 


(13) 
L, + Lins 2 Lin pa. 
L = pls L, L, 5 oe 
2 Lis a 5 Lins L, a Lis 

(14) 


The flux linkage and voltage equations in (1) and (2), in terms of variables referred to the 


stator windings, can now be expressed as 
Raves L, L,, L shes 
Ass (Li i L ie 


Vanes | = Ly, [es + P| aves | 
LV per | = [7 ae a3 P| Avver | 


(15) 


(16) 


Above, the physical voltage equations were given for the induction machine; however, in 
an effort to reduce the complexity of their differential equations we will use a change of 


variables transformation to a new reference frame as described in equation (17). 
ae = K; ca) 
Ree ] = , Fas 
(17) 


where the matrix K, for stator variables and K, for rotor variables is expressed as 


cosO cos(@—-22) cos(@+22) 


K -= sind sin(@-22) sin(@+=) 


boon 
2 2 2 
(18) 
cos cos(B-*2) cos(f+24) 
K,=5 sinf sin(f—2£) sin( +z) 
14 1 
? 2 2 
p=0-6, 
(19) 


In equations (18) and (19), @ is the variable used to transform the f.,. variables for 

which equation (17) are referred. The place holder x in equation (17) represents this 

particular frame of reference that variables are transformed. The place holder f 

represents the voltage, current, or flux linkages in the stator or rotor windings. If no place 

holder is used then the reference frame is considered to be an arbitrary frame of 

reference. When transforming between reference frames equation (20) is utilized where 
y represents the new frame of reference. 

cos(@,-0,) — sin(, - 8, ) 

“K* = -sin(0, - 9, ) cos(0, -0 ) 

0 0 
(20) 


After performing the transformation of equation (17) on each component in the matrices 
on the right hand side of equation (16), the voltage equations can now be expressed in the 


arbitrary reference frame as 


V,, =Vl,, + OA, + PA,, 
=1 1. = OA, + ”,, 
=Nly, + PA, 
vi, =r, +(@-@, )A), + pay, 
vie = ri, —(@-@, Jj, + Ply 
=riy. + PAy, 
(21) 


In equation (21), @.1s electrical angular frequency of the rotor, and @ is the angular 


frequency of the reference frame that the variables are transformed in radians per second. 


Since yw = Aa, and @, = 2760"... , it is convenient to express the inductance component 


of the flux linkages in equation (21) as a reactance instead of as an inductance, beginning 


with equation (22). 
fo) 
Lae = Flas ap Was a Was 
b , 
fo) 
Vas r Las aon car Was 
f , 


o-©o,) 
popes ( r ’ pP ’ 
vag = Pi Loy =F War VE 
b b 
o-o,) 
eg het ( r U p ’ 
Var Play a ee AP —_ Vx 
DM, - 
pone p ’ 
Vo. = "1, g3 —_ Wo, 
@ 


(22) 


The flux linkage component of equation (21) now has units of flux linkage per second 


and represented by the symboly . For ease of understanding the flux linkages in (21) and 


flux linkage per second in (22) are compared in equations (23) and (24). 


Ag = Lig 


lsqs + Lins (i, ole i) 


. . os 
Kas = Lighas 23 Le (i, 7 L,,) 


Ans ? Lilos 
a . . of 
hi = Ligh +L, (i,, + i',) 
pr . . of 
hip i: Lily 27 Las io Tr i.) 
yO f 
Ay 77 Lilo, 
(23) 
. ; of 
Was + Xislas a x ms (i, as L.) 
‘ . of 
Was = X jglas zs X ns (i, os i) 
Wo. ss X los 
po fs . f 
Wir = ibs - a (i,, a i',) 
i fs . 7 
Wy _ Xrly a ae (ix a, ve) 
po poof 
Wo, _ X lo, 
(24) 


The above equations mathematically represent the behavior of a DFIG, and are used in 
the simulation process to arrive at an appropriate physical model that will be discussed 
below. More can be read concerning induction machine theory as used in this thesis 


application in [4]. 


B. SYSTEM CONFIGURATION 


The combined system 1s shown in Figure |. The DFIG is a 175W, 120V, 60Hz, 4- 
pole machine with parameters listed in the Appendix data sheet. The system parameters 
were used in order to effectively model the system in Simulink. The Lab-Volt model 
8231 was used as the DFIG. The stator consists of 3-phase wye-connected windings with 
a turns ratio of Ns/Nr = 10 and resistance R; = 12Q. The rotor consists of 3-phase wye- 
connected windings with resistance R, = 4Q. The reactance of the DFIG are X,,= 180Q 
and X/,=X),= 9Q. 


The DC motor used to emulate the wind turbine 1s the Lab-Volt 8211. It is a 
175W, 120V, 4-pole machine with a nominal operating speed of 1800 RPM. The 
operating speed range used was 1350 RPM to 2150 RPM. 
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The power supply used for the experiments was the Lab-Volt 8821-10. It is a 3- 
phase 120V to 208V AC power supply with a separate VARIAC power supply and 120V 
DC power supply. In order to properly provide power to the Supply-side circuit, the 
VARIAC output is used for the Supply-side SVM with a nominal setting of 60V rms. 
This allows sufficient voltage to be boosted to maintain the DC Link at 200V without 
excessive current transients. Further, a set of 640uH inductors were placed in series (see 
Figure 2.) to the input of the SVM in order to decouple the current sensors in the SDC 
from the switching action of the SVM. 





Figure 2. Simplified Single Phase Circuit Showing the Inductor Connections. 


The mathematical analysis used to arrive at the equivalent simulation matrices 
used in Simulink are shown in equations (25)-(31). The single phase neutral point 


voltages are 


Vn TY ne tM ng = 3, —3N, = 3, volts 


an2 Cc 


(25) 
and the KVL analysis results in 


y 


an2 


+ SL pi, + Rpt, =Voany + (Mh, — My) 
(26) 
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Vino + sb pl, + Kyi, =V 55 (ny) 


(27) 
Von? + sb i, + Kyl, a V sent a3 (n, — Ny) 
(28) 
Where n; and n2 are neutral point voltages. 
Subtracting (26) and (27) gives 
Yip a Vino +sL (i, —I,) +R, —i,) = Ye —V sont 
(29) 
and subtracting (27) and (28) gives 
Vino 7 Ven? + sb (i, —1,) + R(i, =.) = Vi 7 V sent 
(30) 
Defining V4.) Ving =Vae av >Vinz Ven =Vae pes Ad 1, +1, +1, =O the matrix yields 
V 1 -l|li 1 -l|li V 
dc _ab 49 Ly : ae R, : = sab 
Y 4s be 1 2 |i, 1 2 ii, ae 
(31) 


Next, equation (31) can be put into a state space form by defining 
1 -l 1 -l l Vi. an a. V 
N =L, M =R, P= dc _ab sab _ ab 
I. “2 1 2 | Vic be is Vie 
So that (31) can be written as 
Mx+sNx+Pu=0. 


(32) 
Then the state space model becomes 
sx = Ax+ Bu 
(33) 
with 
A=-N'M 
B=-N™'P 
(34) 
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Therefore, the matrices used in Simulink are 
L, L, —L, ' RK, Ky |li, | |b, Ly Ena 
S = — 
i Ly. 21, R, 2R, |/i, L, 2L, 1 || V,. 


The above process represents the common approach used to parallel circuit 


(35) 


analysis and Simulink modeling. Since Simulink is a mathematical based simulation 
program, vice any given circuit solving based program, it is imperative to accurately 
analyze the mathematical characteristics of any given subsystem in order to effectively 


implement Simulink. 


Once an in depth understanding of induction machine theory and system 
configuration was established, specific control methods were explored. The control 
methods included control components, Supply-side control, DFIG-side control, SVM 


scheme concepts, and the encoder implementation. 
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Hl. CONTROLLER SCHEME 


A. CONTROLLER COMPONENTS 


The back-to-back SVM components used for this research are a combined 
package manufactured by SEMIKRON. The package is called the SEMISTACK-IGBT 
and consists of a three phase inverter that 1s coupled to the DC Link capacitor bank. The 
inverters inside the SEMISTACK are made of SKM 50 GB 123D IGBTs that are 
controlled by SEMIKRON SKHI 22 gate drivers. The Supply-side SEMISTACK 1s used 
to maintain the DC Link at 200V and can also be used to control the power factor to the 
grid when in inverter mode. Controlling the grid power factor was not within the scope of 
this thesis, outside of experimentally showing its explicit control, but is a feature that can 
be used to minimize transmission losses in a large scale project. The control for each 
IGBT in the SEMISTACK is implemented through the SDC. Two SDCs are utilized 
(Supply-side and DFIG-side) and include a Xilinx FPGA and a Control board that has 
various inputs and outputs that interface with the Supply-side and DFIG-side circuits. 
One of the inputs 1s an 8-channel, 12-bit A/D converter used to measure the voltage and 
current signals from their associated systems. The DFIG-side circuit also has a 4-pin 
input that samples the rotor encoder pulse waveforms and converts them into a rotor 
position and rotor speed input via an algorithm implemented in SIMULINK. The FPGA 
is programmed to take the sampled inputs and generate the desired gate signals to the 
control board. The control board has BNC (Bayonette Neil-Concealman) connecters that 
connect to the SEMISTACK’s gate drivers. The gate drive signals activate the IGBTs in 
a fashion that will produce the desired 3-phase current in the rotor windings for the 
DFIG-side, or the commanded DC Link voltage for the Supply-side. The procedure for 
using the SDC, architecture of internal components, and the benefits of using FPGA for 


digital control of power systems is discussed in [6]. 


B. SUPPLY-SIDE CONTROLLER 


The Supply-side circuit is powered via a 3-phase AC power supply stepped down 


by a VARIAC to 60V rms. Of note, it 1s important to calculate the DC Link injection 
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voltage as described in [2]. The injected voltage should be less than Vpo ina / 2/2 , 


which corresponds to an injected voltage of 60V rms. There is a set of three inductors 
connected in series to the SEMISTACK, downstream of the Supply-side SDC, which 
decouples the SDC sensors from the switching action of the IGBTs. Figure |. shows, and 
equations (25) through (35) describe their matrix analysis. The SDC senses the stepped 
down power supply line-to-line voltages, Vag and Vgc, and derives the Supply-side 


voltage angle, 0,, as shown in Figure 3. 
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Simulink Model of Supply-side Voltage Angle Calculation. 


Figure 3. 
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Since the SDC only has four voltage inputs it 1s important to be efficient about the 
voltage input configuration. Further, it is not always convenient to input single phase 
values due to the difficulty, at times, to tap a neutral point. This is why it is suitable to use 
the line-to-line voltage parameters as inputs. The derivation of the Vag and Vgc based 


voltage angle calculation as described in [7] and [8] 1s 


Vip Vo ape = Up Ve -(-», = 


(36) 
Note also that i,+7,+i.=0 for a wye connected winding set. 


Placing into matrix form 


(37) 


We define the matrix in dqO terms using 


Vv, =%ly, cos(@)+y, cos(- 27/4) +v, cos(9+27/) 
= ly, cos(@)+y, cos(0-27/\+(-, -v, )eos(0+27/) 
= %)», (cos(#)—cos(0+27/))-+v, (cos(0-274)-cos(@+274)) 


(38) 


and 


Vv, = ly, sin()+v, sin(0-27/) +. sin(0+27/) | 
=Aly, sin(O)+v, sin(0-27/)+(—y, -v,)sin(0+27/) | 
= %)», (sin(@)-sin(0+27/))-+v, (sin(9-24/) -sin(a+2#/))] 


(39) 


After combining terms the line-to-line voltage equations become 


18 


bse) terval 


(40) 
Finally, 0, is calculated using the ArcTan function with v, and vy as the inputs for the 
angle calculation. The ArcTan function gives an output from -27z to 27. Since Xilinx is 
a digital bit-wise based program, it is necessary to implement the last block of Figure 3, 
which converts the ArcTan output to a wrapped 0 to 2’” bit output. The thetaconv2 block 
code is listed in the Appendix. 


The SDC sends the sensed power supply voltages and currents and processes the 
data in order to send the correct gate signal sequence to the SEMISTACK BNC 
connections. The IGBTs will gate on and off and send the supplied voltage to an internal 
capacitor bank, which supplies the DC Link terminals. Of note, each SEMISTACK has 
it’s own capacitor bank, and their associated external terminals, which is used to 
physically couple the Supply-side and DFIG-side DC link to form a common Vpc. The 
last Supply-side circuit input is Vpc , which is sent to the SDC as well. This, of course, 1s 
used to sense and maintain DC Link voltage at 200V. It is also used in a protection 
feature that shuts off the gate signals to the Supply-side SEMISTACK in the case that DC 
Link voltages reach 240V. The implementation of this protection feature utilizes an SR 


Flip-Flop. 


In order to effectively control the DC Link voltage an inner and outer PI 
controller scheme is used. The outer loop controller senses Vpc, as compared to the Vpc 
reference of 200V, and generates a DC reference current (i, ;¢ that 1s sent to the inner 
loop. The inner loop is controlled in the qdO synchronous reference frame so as to easily 
compare i, ,.to the measured i, in DC terms. The derived 6, is now used as the B input, 
equation (19), for the qd0 equation to transform the supplied currents into the 


synchronous reference frame. Figure 4 shows the inner and outer loop Simulink model. 
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Figure 4. Inner and Outer PI Control Loops for the Supply-side Controller. 


C. DFIG-SIDE CONTROLLER 


The DFIG rotor is powered via the DFIG-side SEMISTACK IGBTs through the 
DC Link. The DFIG stator is powered via a 3-phase AC, 120V, 60Hz power supply 
corresponding with a grid. There is a set of three inductors connected in series to the 
SEMISTACK, downstream of the DFIG-side SDC, which decouples the SDC sensors 
from the switching action of the IGBTs. The SDC senses the stator line-to-line voltages, 
Vap and Vgc, and transforms them into the qd0 rotor reference frame using @, (from the 
encoder) as B in equation (19). With the stator line-to-line voltages transformed to the 
dqO rotor reference frame, the DFIG slip angle and frequency can be derived using the 
same method as described 1n the Supply-side controller section using the ArcTan function 
with v, and vg as the inputs for the angle calculation. The DFIG slip angle (6,;;,) 1s 
important since it will be used to control the DFIG in a synchronously rotating dq0 
reference frame, with the g-axis oriented along the stator-flux vector position. This, in 
turn, enables decoupled control of the electrical torque (i,,) and rotor excitation current 


(ig). Traditionally, Os) 1s obtained using the ArcTan function with a rotor referenced Was 
and Was as inputs, along with a offset, in place of the ArcTan calculation using v, and 


vg as the inputs for the angle calculation. The offset is associated with the relationship 
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between current and flux linkage as it relates to voltage. Since the 0,);, calculation used in 
this thesis is based on voltage, no offset is needed. The voltage based calculation of Osjip 
requires a less complex algorithm, and thus less processing, than the flux linkage based 
Os1ip calculation. Both experimental and simulated tests were performed to compare the 
two @.;;, methods. In either test there was no detectable difference in the results for this 
model’s application. Further tests should be performed to determine the extent of 
constraint between the two methods and either’s benefits. Figure 5 shows an overview 
model of the slip angle calculation. The internals of this model can be sufficiently 


described by reviewing Figure 3. 
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Figure 5. Simulink Model of DFIG-side Slip Angle Calculation. 


In addition to the SDC sensing the stator line-to-line voltages, it also senses the 
phase A and B rotor currents (phase C is derived using phases A and B). The rotor 
currents are utilized, along with rotor angular frequency (w,.), to form the current control 
loops for the DFIG. The DFIG PI control loops consist of an outer loop and an inner 
loop. The outer loops takes w, as an input and compares it to a commanded w, ;er . The 
resulting difference passes through the speed PI controller and sends a reference current 
(ig ref) to the current control loop, or inner loop. The reference current passes through the 


current PI controller and compares it to the measured i,,. The resulting waveform is 
pA 


converted to a per unit value and sent to the SVM block, which generates the gate signal 
sequence for the DFIG-side SEMISTACK. The current control loop is used to command 
a specified w, as controlled via iz, or rotor excitation current via ig,. Rotor angular 
position and speed are derived from the encoder, which will be discussed later. Figure 6 


shows the DFIG controller topology. 
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Figure 6. SIMULINK DFIG Controller Topology. 


Both the Supply-side and DFIG-side controllers have unique differences up to their 
respective outputs, however, both controller outputs go to the exact same SVM 


architecture in order to control their commanded parameters. 


D. SPACE VECTOR MODULATION 


The modulation scheme used to create the desired 3-phase waveforms for the 
Supply-side and DFIG-side circuits is SVM. SVM 1s a specific form of Pulse Width 
Modulation (PWM) in which an algorithm involving space vectors are used to control the 
on and off times of pulsed signals. The generated signals then drive the IGBT gate signals 
to the SEMISTACK allowing the user to create a waveform of any magnitude and 
frequency desired. The SVM approach utilized in this controller is obtained from [8-9]. 


The SEMISTACK produces a 3-phase output voltage, as supplied from the DC 
Link, that is controlled by turning on and off its six IGBTs. Figure 7 shows both the 
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SEMISTACK-IGBT configuration and SVM Hexagon for which the on and off 


switching states are based on. The gd-axis is overlaid on the SVM Hexagon and 


represents the axis of the stationary reference frame variables v,andv,, which are 


shown entering the SVM block. The magnitude of the vectors v,,and v,, is equal to the 


vector Vand @ is equal to the ArcTan of the two vectors. 
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Figure 7. SEMISTACK-IGBTs and SVM Hexagon [From 8]. 


There are six sectors on the hexagon and eight possible on and off states that are used to 
produce the V’ vector. The vectors V, and V2, when summed, are also equal to the vector 


V’. In Sector I, V; and V2 correspond to the states (p,n,n)and(p, p,n). The pand n 
correspond to the positive or negative IGBT bus signals for the respective phasesV_ ,V, , 
andV,. The zero vectors, which are represented by the states (p, p,p) and(n,n,n), are 
not shown but are represented as vectors into and out of the page. The magnitude of the 
vectors V; and V2 correspond to the amount of time spent on the switching states in each 
sector. The duty cycles 7, and T, are the times spent on each cycle for the vectors V, 
andV,. The total time spent for one switching period is 7’ 


_ 21 a 


V 
a7 


(41) 
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as 2TV 4, 





y 
2 37 
(42) 
The law of sines on any one of the sectors in SVM Hexagon will produce 
wee VY, 
V3 sin (60° - 0) sin (0) 
(43) 


and substituting equations (41) and (42) into (43) 1s used to find the time of the duty 


cycles for the vectors V, andV, 


mae 








T, = T, sin(60° - 0) 
dc 
(44) 
c= i V3 7 sin (0) 
dc 
(45) 
Wo =1, +1, +1y 
(46) 


SVM has the ability to minimize switching loss and harmonic distortion by controlling 
the order in which the states are applied so that only one gate is fired between sectors [8- 
9]. The switching pattern for each state and sector used in this thesis is shown in Table 1. 
The SVM algorithm is implemented inside each SDC where it is oversampled to 
converge on the performance of a true continuous sinusoidal signal. The digital 
implementation of the SVM scheme used for this thesis is discussed in detail in [8-9] and 
the Simulink diagrams can be found in Appendix C. Figure 8 is a block diagram of the 


functions for each process in the SVM block. 
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SVM Digital Implementation Diagram [From 8]. 


Figure 8. 
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E. ENCODER IMPLEMENTATION 


The encoder used for the DFIG rotor position and speed input is the Microtech 
Laboratory Inc. (MTL) MES-20-200P C. Its implementation is discussed in [10], 
however several modifications were made to commission the encoder in this thesis. The 
encoder uses a single shaft, 200 pulse, open collector output design to produce an Ax, Bx 
and Z,x pulse sequence as seen in Appendix A and Figure 9 (experimental Chipscope 


capture). 


Bo clean state (0/1) 





Tine (3) 


Figure 9. Experimental: Encoder A, B and Z Pulse Configuration. 


Of note, the imperfections shown in Figure 9, and some subsequent figures, are 
caused by a Chipscope plot function artifact due to Chipscope timing discrepancies. 


Pulses A, and By, are in quadrature, and the Z, pulse occurs once every revolution (200 
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pulses) in alignment with the B pulse period. In order to derive the rotor position from 
this encoder input a Simulink block was created to use the A, and Bx, pulses as an 
incremental counter and the Z, pulse as a one revolution reset. Figure 10 is the Simulink 


block implemented for this function. 


Ae 1 
gr 
jj, Paige 





Figure 10. = Architecture Used to Derive Rotor Angular Position from the Encoder. 


The concatenation block takes A,, Ax.; and B, and forms a three bit word and 
relates it to a two bit word (Px and Nx combined). The two bit word is compared to the 
predetermined encoder vector in the ROM block (Px, and Nx is also used to enable the 
counter). The predetermined encoder vector is specific to the encoder characteristics and 
is used to indicate whether the encoder shaft is turning clockwise (binary | 0) or counter- 
clockwise (binary | 1), hence count up or down. A rising and falling edge architecture is 
chosen for the encoder algorithm, which effectively doubles the encoder resolution from 
200 triggers per revolution to 400 triggers per revolution. The two bit word is then passed 
through the Simulink Slice block, which separates the word into two one bit words. Px 
and N; trigger the counter to count up or down with a range of 0 to 2" bits with a wrap 
feature triggered by the reset pusle. Table 2 shows a detailed account of the encoder 


algorithm. 


Zi 


P, Nx | Enable Up | Counter State Encoder Vector 
Logical OR 
ee 





Table 2. Encoder Truth Table with Directed Actions. 


The Z, pulse provides inputs to an AND gate as seen in Figure 10. Since it occurs only 
once per revolution, it is straightforward to see why it is used for the counter reset. The 
counter output is a true account of the encoder position, however the actual rotor 
electrical position is not accurate until the encoder zero point is calibrated to be aligned 


with the rotor electrical zero point. 


The encoder calibration procedure began by a detailed study of the rotors physical 
winding layout. The DFIG is a 4-pole, 3-phase induction machine with 24 slots and no 
accessible rotor windings neutral point. A simple method of aligning the rotor to a zero 
point electrical position would be to supply power to the phase A windings of the rotor 
and stator and allow the magnetizing flux to align the two. Without access to the rotor 
winding’s neutral point, this method became moot. The second, and possibly more 
accurate, method to align the rotor to its zero point electrical position was to supply 
power to the rotor and stator phases via line-to-line. After several experiments with the 
line-to-line method, there were multiple convergences to a rotor zero point electrical 
position to be 318 triggers out of 400. This offset was added to the encoder algorithm as 
the calibration point. Finally, the calibrated encoder signal was passed through a series of 
conversions so that it corresponded to the rotor electrical angular position and also had a 
wrapping range of 0-2'° bits per revolution. Figure 11 shows the experimental results of 
the calibrated rotor electrical angular position. 
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Figure 11. | Experimental: Calibrated Rotor Electrical Angular Position. 


The last input the encoder provided was the electrical angular frequency, w,. This 
derivation was simply a differentiation of rotor electrical angular position with respect to 
time. Due to the limitations of Simulink’s iterative features, the output of this 
differentiation was crude and noisy. However, after passing this output through a first 
order low pass filter with a corner frequency of 10Hz, the signal was surprisingly 


accurate and proved to be very functional for this application. 
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IV. RESULTS 


A. OVERVIEW 


In this thesis, we emulate a Wind Turbine Generator by driving a DFIG via a DC 
motor with variable input torque capability. As discussed in [1-4], we can use a DFIG 
with back-to-back SVMs to accomplish a decoupled Supply-side, and DFIG-side, control 


scheme while allowing power flow to occur between the two systems. 


The two circuits of concern are the DFIG and Supply-side circuits. The DFIG and 
Supply-side circuits are electrically coupled with back-to-back (DFIG-side and Supply- 
side) SVMs that are coupled through a DC Link consisting of a capacitor bank. The back- 
to-back SVMs, with DC Link, provide bi-directional power flow between the DFIG rotor 
and Supply-side power supply. Bi-directional power flow is achieved when the DFIG 1s 
controlled at supersynchronous speeds, and input torque is enough to overcome losses, 


and subsynchronous speeds (Supply-side provides power to rotor). 


B. SUPPLY-SIDE EXPERIMENTS 


The Supply-side power supply (3-phase AC) provides power to the Supply-side 
circuit (stepped down by a VARIAC) and the Stator of the DFIG. The Supply-side circuit 
senses the supply side voltage, current and DC Link voltage and sends the sensed inputs 
to a FPGA. The FPGA is programmed via a Simulation software package called Simulink 
to achieve the desired control of the specified parameters. For the Supply-side circuit we 
use the mentioned sensed inputs to derive control of Vpc so that it is maintained at 200V 
via SVM. There is also a built-in protection feature in the Supply-side circuit that will 
turn off the SVM if a VDC of 240V is reached. Within the Supply-side controller 
architecture, however, there are parameters that can be used to control the system power 
factor. The research of this feature, and its effects on the grid, was beyond the scope of 
this thesis, but some simulation and experimental tests were conducted to show the 
performance of the controller. If a system, connected to the grid, has explicit control of its 


power factor, it can be used to reduce transmission line losses by controlling the amount 
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of reactive power being distributed. Figures 12 and 13 show the simulation and 
experimental steady state shifts in power factor when ig ;er was commanded at 2A and - 


2A. 
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Figure 12. — Simulation: Supply-side Voltage and Current with ig ;er at 2A and -2A. 
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Figure 13. 
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Experimental: Supply-side Voltage and Current with ig ,.¢at 2A and -2A. 


C. DFIG-SIDE EXPERIMENTS 


The DFIG circuit power supply is provided by the DC Link and DFIG-side SVM. 


The DFIG circuit senses rotor current, stator voltage, rotor speed, and rotor electrical 


position (as provided by an encoder). These inputs, like the Supply-side circuit, are then 
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sent to an FPGA, which is programmed to control the DFIG parameters. The DFIG 
circuit will take the rotor speed input and derive a current reference that is used to control 
the rotor speed using a PI controller. It 1s important to show that we can control the DFIG 
rotor speed since traditional Wind Turbine Generators would use turbine blade pitch 
control in tandem with a synchronous generator to provide a constant speed, and hence a 
constant output frequency. Further, it 1s noteworthy that a wind-speed profile can be 
adapted to a specific Wind Turbine (outfitted with a DFIG) in order to program the DFIG 
to run at specified speeds for a given wind speed. This tool can be used to provide an 
optimized rotor speed for a given wind speed. For our application we will show that the 
DFIG rotor speed can be controlled (subsynchronous, synchronous or supersynchronous) 
for any given input torque and thus distributes the power generated by the DFIG as 
predicted between the rotor and stator. Further, this experiment will demonstrate the bi- 


directional power flow of the back-to-back SVMs showing power recovery via the DFIG 
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Figure 14. = Simplified Circuit of Wind Turbine DFIG System. 
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As illustrated in Figure 14, the Supply-side SVM converts a 3-phase AC source to 
a constant DC Link voltage of 200V. The DC Link voltage is then inverted into a 3-phase 
waveform via the DFIG-Side SVM (during modes for which power is being delivered to 
the DFIG). The digital control for the Insulated Gate Bipolar Transistor (IGBT) network 
is implemented via the FPGA. The FPGA interface samples the rotor current and the 
stator voltage of the DFIG through an A/D converter. There is an input from the encoder, 
0,, to provide the position and speed of the rotor to the FPGA through the interface circuit 
board. An algorithm then calculates the slip frequency of the DFIG, which is used to 
derive the needed rotor current frequency and amplitude to maintain a specified rotor 
speed, and therefore a controllable output frequency. Since the Supply-side and DFIG- 
side controller topologies are nearly identical, it follows that the DFIG-side controller can 
control system power factor in the same way that the Supply-side did. Figures 15 and 16 
show the simulation and experimental steady state shifts in power factor when ig ;er was 


commanded at 2A and -2A. 
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Figure 15. Simulation: DFIG-side Voltage and Current with ig ;e¢at 2A and -2A. 
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Figure 16. | Experimental: DFIG-side Voltage and Current with ig ,.¢at 2A and -2A. 


Figure 17 demonstrates an input torque transient representing a large wind gust. 


By commanding a specified supersynchronous speed, while taking into account windage 


and friction losses of the system, it 1s clear that the input torque transient produces a 
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power flow reversal in the rotor. In other words, the rotor mode transitions from power 


flow from the rotor (negative rotor current) to power flow to the rotor (positive rotor 






F : ie 


Figures 18 and 19 (experimental results) parallels this transition with the steady 
state captures corresponding to the simulation transient. The experimental results were 
run through a low pass filter with a corner frequency of 80 Hz in order to create a useful 
plot of the event. Due to windage and friction losses, a commanded supersynchronous 
speed was ordered to achieve the power flow transition results during the constant rotor 


speed experiment. 
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Figure 18. | Experimental: Rotor Current and Voltage with Power Flow from the Rotor. 
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Figure 19. — Experimental: Rotor Current and Voltage with Power Flow to the Rotor. 


Additionally, Figure 17 (simulation model) shows a smooth transition as the 
power flow reverses between power flow from the rotor to power flow to the rotor. This 
transition is modeled using a constant commanded rotor speed along with a 10% decrease 


in input torque. 


Figure 20 demonstrates a large, constant input torque representing a constant wind 
speed. By increasing the commanded rotor speed by 30%, while taking into account 
windage and friction losses of the system, it is clear that the rotor speed transient 
produces a power flow reversal in the rotor that corresponds to subsynchronous, 
synchronous and supersynchronous speeds. In other words, the rotor mode transitions 
from power flow to the rotor (positive rotor current) to power flow from the rotor 


(negative rotor current). 
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Figure 20. | Simulation: Increase in Rotor Speed with a Constant Input Torque. 


Figures 21 through 23 (experimental results) parallels Figure 20 with the steady 
state captures corresponding to the simulation transient. The experimental results were 
run through a low pass filter with a corner frequency of 80 Hz in order to create a useful 


plot of the event. 
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Figure 21. | Experimental: Rotor Current and Voltage While Subsynchronous. 
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Figure 22. 
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Experimental: Rotor Current and Voltage While Synchronous. 
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Figure 23. | Experimental: Rotor Current and Voltage While Supersynchronous. 


In conclusion, we showed that we can control the DFIG rotor excitation to 
provide a constant output frequency. Our application demonstrated that the DFIG rotor 
speed can be controlled (subsynchronous, synchronous or supersynchronous) for a given 
input torque and thus distributes the power generated by the DFIG as predicted between 
the rotor and stator. Finally, this experiment also demonstrates the bi-directional power 


flow of the back-to-back SVMs showing power recovery via the DFIG rotor. 
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V. CONCLUSIONS AND FUTURE RESEARCH 


A. CONCLUSIONS 


The Wind Turbine Power Generation Emulation via Doubly Fed Induction 
Generator Control System was designed, simulated, constructed and tested using standard 
engineering principles and theory. Since the research and simulation phases were aptly 
utilized, the development and construction phases were relatively inexpensive and short. 
A series of test were conducted and validated the design. The experimental results 


showed that the system responded according to the simulation. 


With a baseline prototype that functions correctly, the Naval Postgraduate School 
now has a means to build a large scale prototype for laboratory purposes and future 
research. The future research on this topic will be above and beyond current publications, 


and will have an enabled freedom to exploit the abilities of the DFIG. 


B. FUTURE RESEARCH 


A third SDC should be used to take control of the DC motor. This Wind-side SDC 
can be programmed to mirror current wind speed profiles that are available from the 
National Renewable Energy Laboratory. With a wind speed profile controlling the input 
torque of the DFIG, the DFIG-side SDC can be slightly modified to include a rotor speed 
lookup table. This lookup table can be used to derive an optimal rotor speed for a given 


wind speed, and thus increase the efficiency of DFIG power extraction. 


Second, DFIG-side feed forward terms, based on slip speed and flux linkages, can 
be implemented as described in [1]. The feed forward terms should lower the DFIG-side 


PI controller time constant without any detrimental effects on stability. 


Finally, Supply-side feed forward terms, based on the stators synchronous angular 
frequency and stator winding inductance, can be added. Again, these terms should lower 


Supply-side PI controller time constants without any detrimental effects on stability. 
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APPENDIX A: DATASHEETS 


SEMISTACK - IGBT 


SEMITRANS Stack”) 





Three-phase rectifier + 
inverter with brake 
chopper 

SEMITEACH - IGBT 

SKM 50 GB 123D 

SKD 51 

P3/250F 


anes == womens pe 


* Mult-function IGBT converter 


General dimensions 


This technical information specifies semiconductor devices but promises no 
actenstics. eae ae ga alla made regarding 
ertorn Suitability. 


BéU + B6C! + E1CIKF 


Power Electronic Systems — SEMISTACK 08-06-2005 @by SEMIKRON 
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Model 8231 Three-Phase Wound-Rotor Induction 
Motor 





THA PRESS -SOn oe Pe a ee 








Each phase of the stator windings of this motor is inde- 
pendently terminated and identified on the faceplate to 
pennit operation in either delta or star (wye) configura- 
tion. The rotor windings are brought out to the faceplate 
via extemal slip rings and brushes. This machine can be 
used as a wound-rotor induction motor, phase shifter, 
single-phase vanable coupling transformer, three-phase 
transformer, selsyn control, frequency converter or 
asynchronous induction generator. The speed of this 
machine can be controlled through the use ofthe Three- 
Phase Rheostat (Model 8731). 


SPECIFICATIONS 
Model 8211 DC MotoriGenerator 


Motor Output Power] 175 W 


14.1 at ib) 


tion Motor JAN380 V - 50 Hz 240/445 V - 50 Hz 
femrtegeees fey [ey ey 


Dimensions (H = Wx D) | 308 x 2 4 Tam { 121 x ix 173 in) 
Net Weight) 13.5 — (207 Ib} 
3 waaneo V-s0z | 24015 V— 90 i 


— 

peed] —1800rmin | 280 rimin | STS emi 

Ful Load Curent 
Dimensions (H x Wx D) |) 308 x 287 x 440 mm (12.1% 11.5% 17.3 inj 


Net Weight 
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POWER SUPPLY 


AATIPG 








a 0) Tie Hz 
POWER SUPPLY 
MODEL 8821 
The Power Supply provides fixed and vanable AC A voltmeter, connected through a selector switch, 


and DC voltage sources, all terminated by color-coded monitors the variable AC and DC outputs and fixed DC 
4mm safety sockets. Independent circuit breakers, reset output. A 24 V AC output provides a low-voltage supply 
atthe front panel, protect the input to and output from the required to operate other EMS equipment such as meter- 
Power Supply. Indicator lamps monitor the presence of ing modules and modules used in the Power Electronics 
input voltage in each phase. When a phase leg of the Training System. 

site's power service is out, the lamp goes off to reflect this 

condition. 





SPECIFICATIONS 


Model £621 —- Power Supply 120/208 V— 60 Hz 220/330 V — 30 Hz 240/413 V — 30 Hz 
! 120/206 V 220/360 ¥ 24415 V 


service installation| 20 A, aA Saas 5 wines, star (Wiye} connected, including neutral and 


Fimed AC 3-Phase == 15 A S80 V— 104 240/415 ¥—- 104 
Vanable AC 3-Phase| 0-120/208V-354 0-220/380 V-3 A A404 1SV-3A 


Vanable DC 0-120 V-38 4A 0-220 V-35A 0-240 V-3.A 
Fmed DC 120 V-2A 220V-14 240 V-1A 


Low PowerAC|] 24V-34 


Wall Gutlet included NEMA L21-20 NEMA L22-20 CLIPSAL 5650520 
3m (10 ft 


Dimensions (Hx Wx D)| 308 x 287 x 500 mm (12.1% 11-3 x 18.7 in) 


Net Weight) 13.4 kg (40.5 Ib) 
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A, B, 2 phase 
A, B phase SIN wave, 2 phase square wave 
SIN wave 1.5 Vp-p0.9 V offset 2.0VE0-2V 





A, B, Z phase output Onarp oulput current SrA Max. 
Harmonic -<istortion factor to be within 10% 
(Mewsuring cond itiog to be within 20 iz, 
atfecths valea maar distortion factor mane Be instrument 
SOkHz 
Oxon pha tenes | ee acne 
Starting torque 2010 3N-m (20gf-om) of less 
iicumbinicadct| Fimiiel  47NGLSkg 
shaft (etectical | Thrust 4.9N (0.5kgf) 
Working ambient temperature | o°C~50°C 
AK35%~90%¢ no dewing 
Storing ambient tompernture -20°'C~80°G 
iiientinn seciaheinee Durability 55Hz, double amplitude 1.5mm 


2 hours mach in xX, ¥, and 2 directions 
Durability 500m/s? (about 50) 


ten ee 3 times cach in X,Y, and 2 directions — 
ineulated shield cable (lengtty 1m). . 
Mass 70g 
! Output circuit diagram 








A capacitor (0.1 uF} is connected between OV and FG (frame ground). 


| Output waveform 


pasleghirabnipa adept 


Tn ‘Aivtsions 
; UB AST 
gett 


MWS. cer hea 





Arpiaae ears. 
Tee moeition of 2 pire against A, B phase is not specified. 
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| Specifications Built-in muttiplication circuit (X2+4-><8-<16) 

































































Supply volage valet Com oben OCSY “Sia 1086 
__Currentconsumption | GOmAorless (underno load) 
Detection system | increments) 
| Qutpart in plo ean Ex ey )} 
{Pulse number/rotation) Bs00xi6 (40,00 
Output phase A, B, Z phase 
Output capacity | Open calector out putload voltage DG13.2V max 


Line driver outpul!50kH2% (by multiplication) 
Macimnum reaponse BeUENCY) \oHace cutput-Open collector output: 100KHz 







































Starting torque 2x10°oNem Scoala oh less 
Allowable oad ot|  Flaciial 14.7 (1.Skg 
alata 4.8N (0.5kg!) 
eae ee : 
j -10'C~70'C 

dont ___ FRE E-807E no dewing 


=20'°C—s0°C 
Durability 55Hz, double amplitude 1.4mm 
2 hours: each in X, ¥, and 2 directions 


Durability 500mvs? (about SOG) 
| 3 fimes each in X, ¥, and 2 directions 


Outside dlamotor $4.2 Score vinyl wire 
Insulated shield cable (length 1m) 


ibs I 





| nL 
A capi 1004 ail) in connected between OV and PG (ieee qeounat) 


| Output waveform 
CW rotation (CW rotation as seen from ft surface 
T 


ae 


f phase 
|___TaT2 


T = T#0.3T 
i= i/4 T+6.51 
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APPENDIX B: MATLAB M-FILES 


A. MATLAB INITIAL CONDITIONS FILE 


DFIG.m 


Vdc=200; 666DC Link Voltage Setpoint 

V_phase=sqrt (2) *120; SsSupply-side voltage 

Ks=2/3" [1 =172 =17250 =soqrc(o)72 sort (2)-727172 172 17215 Sabc to gqd0d 
transformation in the stationary frame 
f fund = 60; Ssscbase frequency 


o\? 


6B 


o\? 


omega_b = 2*pi*60; 
oversample=2; 


step_ct=1; 

pulsect=1800/step_ct; 

clock_freq=25e6; 6s6eclock frequency 

tstep = 40e-9*step_ct; %%S%step size 

tstop=4; 

666656556665 SSSSSS566Stator PI Gains$$%%%SSS5S5S5S5S5S5SSS55S SSS 
Kp_vdc= 


Ki_i_s_icq=0; 
Ki_i_s_icd=0; 
L£=300e-6; 6%6$choke inductor 
Rf=.05; 
Amat_indI = -inv([Lf -Lf;Lf LfFt+Lf)])*RE* [1 -1;1 2]; 
Bmat_indI = -inv([Lf£f -Lf;L£E LE+LE]); 
Ciat nd = ib O30: LL gs =b 1s 2lc = -la-Ib 
Dmat_indI = zeros(3,2); 
one_zero_state=0; 6Set to one so that only one zero state is used 
in modulation 
if one _zero_ state == 1 
gainl = i; 
gain2 = 0; 
else 
Sain 1/2. 
gainz = 1; 
end 


twopiby3 = 2*pi/3; 
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poles = 4; 
polesby2J=poles/2/(.089) ; 
TL=200/ (2/poles*omega_b) ; SSSDFIG base torque calculation 


SESESSSESSESSSESROtCOK Speed/Current PI Gains%%SSSSSSSSSS% 
Kp_1=80; 

Ki_3=100; 

Kpgain_speed=.0333; 

Kigain_speed=.0033; 

CESSES SSSCCCE EES SS SESC EEE ESE SSS ECE C EEE SSE SE CECE EEE ESSEEESE ESE SS 


rs=12; S66Sstator res 
rr = 4; S66Srotor res 
Xls =9; 
Xm =180; 
Ar = 23 


D= (X1s+Xm) * (Xlr+Xm) -Xm%*2; 
rsbyXls = rs/Xls; 

rrby xr -= rr/ x17; 

Xagq = 1/(1/Xm+1/X1s+1/X1r) ; 
Xad = Xaq; 

XagqbyXls = Xaq/Xl1s; 
XagbyXlr = Xagq/Xlr; 
XadbyXls Xad/Xls; 
XadbyXlr = Xad/Xlr; 
XagqbyXm = Xagq/Xm; 

XadbyXm = Xad/Xm; 
pS1_9qs1c=0; 

DS1.-ds1¢c=0;7 

psi_gric=0; 

psi_dric=0; 

omegar_ic = omega_b; 
Kp_igs=50; 

Ki_igs=.1; 

Kp_igr=50; 

ha or. 

Fomet = (0 0 0 tel: Lt 2 072 23 073 3.0. 0] 


$66S56565655S6S656S56S6S6SCENCOdeEr Speed 
VeEctorsssss sss 6666S SESS SS SSCS EES SS SSESEEE SESE SS 
Index=[1:2%°12]; 

reciprocal=2%-6./Index; 

$$ $SSS$SSSSSSSSSSS$SENCOdEr 

VECtoOrsssss SSS SSS S6SCCC CEES SESE SCEE EEE SSS ESEEE EES ESSE SS 


Feee 


output_vec=[ 


0 
Oe aes 
2;...eMotor is turning in the CCW falling edge 
1;...sMotor is turning in the CW falling edge 
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oP ol? 


oP ol? 
oP ol? 


oP ol? 


“we 


~e 


.sMotor is turning in the CW rising edge 
.sMotor is turning in the CCW rising edge 
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B. MATLAB M-FILE USED FOR SPACE VECTOR MODUALTION 


cverflow3.m 

function [sectorl, sector2, sector3, sector4, sector5, sector6, 
overflow3 (x) 

$gain = xfix({xlUnsigned,10,7},2.359296/3);%for 60 hz 


Gain = xtix({xlUnsigned, 10,7},2.359296); stor leo hz 
stempv=gain*x; 
tCempv=x; 
if tempv<=171-1 
sectorl=xf1ix({xlBoolean},1) 
sector2=xfix({xlBoolean},0); 
sector3=xfix({xlBoolean},0); 
) 
) 


7 


sector4=xfix({xlBoolean},0); 


sector5=xf1ix({xlBoolean},0 
sector6=xfix({xlBoolean},0); 
Z=xfix({xlUnsigned,10,0},tempv) ; 
elseif tempv<=2*171-1 
sectorl=xfix({xlBoolean},0O) 
sector2=xfix({xlBoolean},1); 
sector3=xfix({xlBoolean},0); 
) 
) 


7 


7 


sector4=xfix({xlBoolean},0O); 


sector5=xf1ix({xlBoolean},0 
sector6=xfix({xlBoolean},0); 
zZ=xfix({xlUnsigned,10,0},tempv-171); 
elseif tempv<=3*171-1 
sectorl=xfix({xlBoolean},0O) 
sector2=xfix({xlBoolean},0O); 
sector3=xfix({xlBoolean},1); 
) 
) 


7 


7 


sector4=xfix({xlBoolean},0); 


sector5=xf1ix({xlBoolean},0 
sector6=xfix({xlBoolean},0); 
Z=xfix({xlUnsigned,10,0},tempv-2*171); 
elseif tempv<=4*171-1 
sectorl=xfix({xlBoolean},0O) 
sector2=xfix({xlBoolean},0O); 
sector3=xfix({xlBoolean},0); 
) 
) 


v 
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sector4=xfix({xlBoolean},1); 


sector5=xf1ix({xlBoolean},0 
sector6=xfix({xlBoolean},0); 
zZ=xfix({xlUnsigned,10,0},tempv-—-3%*171); 
elseif tempv<=5*171-1 
sectorl=xf1ix({xlBoolean},0O) 
sector2=xfix({xlBoolean},0O); 
sector3=xf1ix({xlBoolean},0); 
) 
) 


7 


7 


sector4=xfix({xlBoolean},0); 


sector5=xf1ix({xlBoolean},1 

sector6=xfix({xlBoolean},0); 

zZ=xfix({xlUnsigned,10,0},tempv—4%*171); 
else 

sectorl=xfix({xlBoolean},0); 
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Z. 


sector2=xfix({xlBoolean},0); 
sector3=xfix({xlBoolean},0); 
sector4=xfix ({xlBoolean},0O) 
sector5=xfix({xlBoolean},0); 
sector6=xfix({xlBoolean},1); 
zZ=xfix({xlUnsigned,10,0},tempv—-5*171); 


7 


end 


ramp2mod.m 


function z = ramp2 (x) 
gain=xfix({xlSigned, 20,19},1/1800) 
Z=xT 1x ({xlSigned, 14,13)},x*Ggain) ; 


thetaconv2.m 


function [y] = thetaconv (x) 

gainl = xfix({xlSigned,14,10},2*3.14); 

gain2 = xfix({xlSigned,14,10},1/gainl1) 

1f x<0 
y=xii1x({xLUnsigned, 10,0}, (x+gainLl) *gain2*1024); 
else 

y=xfix({xlUnsigned,10,0},x*gain2*1024) ; 

end 
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C. MATLAB M-FILES FOR CHIPSCOPE INTERFACE 
Black_box_dc_machine2_config.m 


FURST code contig(this block) 


o\? 


Revision History: 


o\? 


o\? 


18-Dec-2008 (152715 hoeure) = 
Original code was machine generated by Xilinx's System 


o\? 


Generator 
% after parsing H:\Docs\work_files\Docs\classes\EC3130_2009\DC 
% machine lab\black_box_dc_machine2.vhd 


this_block.setTopLevelLanguage('VHDL'); 


this_block.setEntityName('code'); 


& System Generator has to assume that your. entity has a 
combinational 

% feed through; 

% if it doesn't, then comment out the following line: 


this_block.tagAsCombinational; 


this_block.addSimulinkInport ('ind'); 
this._block.addSimulinkinport ("ale clock’); 
this _block.addSimulinkinport (*indZ*)-; 


this block.addSimulinkOutport (*ouEd" ); 
this_block.addSimulinkOutport ('open_loop'); 
this_block.addSimulinkOutport ('speed'); 


eutd port = this_block.port (*outd"); 
outd_port.setType('UFix_1_0"'); 

open_loop_port = this_block.port ('open_loop'); 
open_loop_port.setType('UFix_1_0'); 

speed_port = this_block.port('speed'); 
speed_port.setType ('UFix_8_0'); 


if (this_block.inputTypesKnown) 
%& do input type checking, dynamic output type and generic setup in 
% this code block. 


at (this_block.port (*ind*)..w1dth ~= 1); 
this_block.setError('Input data type for port "ind" must have 
width=1.'); 
end 


this_block.port('ind').useHDLVector (false); 
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af (this _Dlockwport (*2la cleck*) .width ~= 1); 
this_block.setError('Input data type for port "ila_clock" must 
have width=1.'); 
end 


this_block.port ('ila_clock') .useHDLVector (false); 


zt (this_pbleck.pert (*ind2*) width ~= 46); 
this_block.setError('Input data type for port "“ind2" must have 
width=48.'); 
end 


if (this_block.inputRatesKnown) 
setup_as_single_rate(this_block, 'clk','ce') 


{e) 


end % if (inputRatesKnown) 


&% Add addtional source files as needed. 


I Na | Add files in the order in which they should be 

% compiled. | If two files "a.vhd" and "b.vhd" contain the entities 
& entity_a and entity_b, and entity_a contains a | component of type 
& entity_b, the correct sequence of | addFile() calls would be: | 

— this block.addFPile("*b.vhd"); | this _block.addFile('a.vhd'); 

% this: block.addFale(*")¢ this block.addFale(**)-; 


this _block.addFile('black_box_dc_machine2.wvhd'); 


return; 


function setup_as_single_rate (block, clkname, cename) 
inputRates = block.inputRates; 


uniquelInputRates = unique (inputRates) ; 
if (length (uniqueInputRates) ==1 & uniquelInputRates (1) ==Inf) 
block.setError('The inputs to this block cannot all be constant.'); 
return; 
end 
if (uniquelInputRates (end) == Inf) 
hasConstantiInput = true; 
uniquelInputRates = uniquelInputRates(l:end-1); 
end 


if (length(uniquelInputRates) ~= 1) 
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block.setError('The inputs to this block must run at ae single 
rate.'); 


return; 
end 
theInputRate = uniquelInputRates (1); 
for 1 = 1:block.numSimulinkOutports 
block.outport (i) .setRate (theInputRate) ; 
end 
block.addC1lkCEPair (clkname, cename, theInputRate) ; 
return; 
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APPENDIX C: SIMULINK/ XILINX MODEL OF WIND TURBINE 
GENERATOR SYSTEM (DFIG-SIDE) 
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A toD conversion 
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thetar_out 


Encoder1 
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alt double 
double double 


Counter CONVSTn ai 
Constant Relational 
a ve? 
double » C) 
Constant2 Relational’ 
CMult2 


d double 









—_—__r 
ch3 
CurrentState pouble 
double double 
Mealy State Machine double 
a CMult3 
doub : double 
double ' 
z —— 
ch2_ 4 
DATA IN - 
= data_rate > 
Delay! 
delay added for metastability 
Subsystem 
Constant Gain3 
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AD state 
double 
double 
" double ; 
double ———>(5 ) 
data_rate 
Counter1 : 
Constant4 Relational3 


AID CONVERTER (VOLTAGE) 
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double 


out 
doubld double 
double “ot 5 
Counter CONVST2n double 
Constant! Relational —>(4) 
ch4 


e CMult2 





Delay2 Mealy State Machine 





AD state 










data_rate 
Subsystem 
double 4 | double 
Z 
doiiie DATA _IN2 
Delay1 
delay added for metastability 
double Multiport CMult3 


Switch 
lbs 
Gain! double 
4 
chi 


AID CONVERTER (CURRENT) 
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ch2 


Positive (Pk) and Negative (Nk) bits are created 
based on Ak transitions and Bk 


Idc 






double 















rad/s conv\ 
IIR first order filter at 10 Hz 
Slice Pk 
tl double 
Pk slice 
double matt 7 double E 
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V_ab to V_qd1 


double 


Vab_s 


Vab_s 


double 


Vbc_s 


double 


thetar 


rotor 
controller 









Inverter4 
double 


ROTOR CONTROLLER 
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double 






double 







double 


SineCosine2 






AddSub2 


double : ( , ) double 
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Negate' 
double 
double 
SineCosine1 
SLIP DERIVARTION 


double 






double iy s double 
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double AddSub1 





double double 









CORDIC ATAN 
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double = 


thetay 
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offset AddSub3 
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double | 
wr ref EB double | 
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APPENDIX D: 


SIMULINK/ XILINX MODEL OF WIND TURBINE 
GENERATOR SYSTEM (SUPPLY-SIDE) 
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APPENDIX E: TRANSFORMATIONS 


Vip =Vq — Vp Myo =V, —V. =V, —(-v, —¥, )3 where v, + v, +v, =0 for a wye connected capacitor bank (with the neutral floating) 


Note also that i, +i, +i, =0 for a wye connected transformer winding set. 
pth be beste be] 
= => =_— 
Ve, Le 2h Veli opal: VE 
v, = Aly, cos(0)+v, cos(O-27/) +, c0s(0+2%/) | ref. Krause text 
= Aly, cos(0)+v¥, cos(0-27/')+(—y, —v, Jeos(9+27/) 





3 


“3.3 


-l | -1 1 


xa ee ORE Wek egy ee ETE 


(0,)  sin(@, feed 


Vin Lae 


-cos(9, +7) —cos (0, ) 


; transformation used for measured line-to-line voltages 





in(, ) —cos(8, +0, 
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